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^ Abstract. 

In Hele-Shaw flows at vanishing surface tension, the boundary of a viscous 
OQ fluid develops cusp-like singularities. In recent papers [TJ [5] we have showed 

^ that singularities trigger viscous shocks propagating through the viscous fluid. 

Here we show that the weak solution of the Hele-Shaw problem describing viscous 
(^ shocks is equivalent to a semiclassical approximation of a special real solution 

^^ of the Painleve I equation. We argue that the Painleve I equation provides an 

lO integrable deformation of the Hele-Shaw problem which describes flow passing 

^^ through singularities. In this interpretation shocks appear as Stokes level-lines 

of the Painleve linear problem. 



1. Introduction 

Hele-Shaw flow describes a 2D viscous incompressible fluid with a free boundary 
at low Reynolds numbers. The fluid is driven to a drain by another, inviscid, 
incompressible liquid or a gas (see e.g., [3]). Under this process the boundary of 
the fluid evolves in an unstable manner, revealing singularities occurring at a finite 
time. 

The governing equations - Darcy's law - can be easily derived from Navier-Stokes 
equations under the assumption of a Poisseuille profile. In units such that the fluid 
density and hydraulic conductivity are set to one, velocity of the fluid is proportional 
to the gradient of pressure 



-Vp. 



(1) 



In incompressible fluids, V ■ v = and pressure is a harmonic function, Ap = 0, 
everywhere except at a drain (a marked point set at infinity), where the fluid 
disappears with a constant flux Q = § \ x d£ i.e. the amount of fluid removed 
in the time interval dT is QdT. 

In the regime where the surface tension at the fluid boundary may be neglected, the 
case we consider, pressure is a constant p\r = along the fluid boundary T. 




Figure 1. A finger growing in time (sequence of blue lines) toward the (2,3)-cusp 
singularity. The dashed line is the branch cut of the height function. 



The Hele-Shaw problem is just one particular realization of a broad class of growth 
processes governed by the Darcy law. Darcy law states that the normal velocity of 
a line element of the boundary di is proportional to its harmonic measure density. 



Harmonic measure is the probability of a Brownian excursion starting at a marked 
point (in this case it is a drain) to exit the domain through that boundary element. 
The probability density of Brownian excursion is the Poisson excursion kernel, equal 
to the gradient |Vp| of a harmonic function p with a source at a marked point and 
Dirichlet boundary condition. 

A well-known stochastic realization of the Darcy law is DLA - diffusion-limited 
aggregation [3] . It is realized through Brownian excursions of particles with a small 
size h emanating at a constant rate from a distant point, with absorbing boundary, 
such that the stopping-set cluster grows. The Darcy law emerges at vanishing particle 
size, /i — !■ 0. 

Computer experiments with DLA |3] show that a boundary, initially featureless, very 
quickly develops into a branching graph with a width controlled by the size of one 
particle h. This signals that the limit /i — )■ is impossible. 

Darcy's law ([T| also indicates that Hele-Shaw flow tends to a microscale. Non-linear 
differential equation ([I]) is ill-defined. A smooth initial boundary first evolves into 
a fingering pattern, then fingers at a finite, critical time develop cusp singularities 
[3 El [71 [8] . At that point the differential form of the Darcy law stops making sense. 
This phenomenon constitutes the major problem in the field. 

Unlike in fluid dynamics, DLA processes possess a dimensional parameter - the 
particle size h. It plays the role of a minimal area, which regularizes singularities 
emerging in fluid dynamics. 

Comparing Hele-Shaw flows to DLA suggests that, after a singularity is reached, the 
flow should feature shocks - a graph of curved lines where pressure (and velocity) 
suffer discontinuities. Such solutions of ill-defined non-linear differential equations 
are known as weak solutions [9]. 

In Refrs. [H [2] we developed the weak solution of Hele-Shaw problem, which we 
believe to be a solution of the problem of Hele-Shaw singularities. Within this 
solution, a singularity triggers a branching graph of shocks or cracks propagating 
through the viscous fluid. Shocks at small Reynolds number is a new phenomena, 
not yet being observed experimentally. We refer them as viscous shocks. An emerging 
pattern of viscous shocks is reminiscent DLA patterns. 



The weak solution of the Hele-Shaw flow is based on an integral form of the Darcy 
law. This form of the non-linear equation suggests an interpretation of the Darcy law 
as an evolution of a real Riemann surface (a complex curve with reality conditions). 
A singularity signals that the evolving Riemann surface changes its genus. 

A similar problem is known in a different field - a slow dynamics of modulated 
fast oscillatory non-linear waves (see e.g., [lO]). If a non-linear wave features a 
separation of scales, in the form of a slow modulation of fast oscillations, the slow 
part is described Whitham equations pTl 112] . Conversely, the Whitham equations 
can be interpreted as the evolution of Riemann surface^ 

In Refrs. [H |2] we described hydrodynamics of shocks initiated by the most generic 
cusp singularity, referred to as a (2,3)-cusp. In local Cartesian coordinates aligned 
with a cusp axis, the shape of a critical finger is y^ ~ x^, as in Fig. [l| 

Right before the occurrence of a (2,3)-cusp singularity, a Riemann surface describing 
Hele-Shaw flow is a sphere. Right after a singularity, the sphere transforms 
into a torus (or elliptic complex curve). Exactly the same transition appears in 
a semiclassical approximation of special, often called "physical" solutions of the 
Painleve I equation [T3| fl^. 

Appearance of the Painleve equation in Hele-Shaw dynamics is not accidental. In 
this paper we clarify this relation. We show that the Hele-Shaw singular flow and 
emerging shock pattern are directly related to "physical" solutions of Painleve I 
equation: 

Utu - 12uut = h, (2) 

where the Painleve -function u{t) is a real function for real values of the 
argument. 

As it is seen in Fig. |2j the Painleve -function behaves drastically different at large 
negative times t < and large positive times t > 0. At t < the Painleve -function 
is smooth. It turns to an oscillatory regime at t > 0. 

The Painleve -function may be linked to the position of the tip of the finger e(T). 
The relation changes while passing through the singularity: at the smooth phase 

I An observation that a Hele-Shaw flow has a form of Whitham equations appeared in Ref. [T5] . 




Figure 2. Illustration of a physical Painleve -function by means of Whitham 
approximation: at t < the Painleve -function is smooth u sa ^J —htj^. It 
features fast oscillations at i > 0, but the change of the period and amplitude 
remain smooth. The (orange) dashed line is the envelop of minima u(i„) ~ 
0.553594 Ay/it„/6. The picture captures only the asymptotes at large |t|, not a 
crossover between smooth and oscillatory regimes. 



(large negative T) the tip progresses as e(T) = — 2n(t) < 0, T — )► — oo; at large 
positive T (after the singularity) the tip retreats, its position follows a slowly moving 
position of minima of the Painleve function e(T) = ui^^) < 0, where t„ is a sequence 
of time, labeled by an integer n, such that u{tn) = 0. In both cases, the relation 
between a "slow" time T of the Hele-Shaw flow and a "fast" time t - an argument of 
the Painleve function is 

T=-, T<0; T=-^, T>0. 3 

o o 

The Hele-Shaw flow is directly related to the semi-classical {h — )■ 0) limit of 
the Painleve "wave-function" \l/ (X, t) - solution of the linear ordinary differential 
equation associated with Painleve I equation: 



d^ - 2u{t) 



^ = X^. 



(4) 



The argument X = x + iy (complex spectral parameter of the linear problem) plays 
a role of a complex coordinate of the fluid. Shocks appear as Stokes level lines of the 
Painleve wave-function. Hydrodynamic objects are obtained from the wave-function 
by averaging over fast time oscillations. In particular, pressure reads 



p{x,y) = — limRe /i5rlog\E'(x -|- iy,t) = —6 limRe dt\og'^{x + iy,tn), 

where overline means an average over fast time oscillations. Alternatively one can 
evaluate the wave-function at discrete time moments (tsl) . 



We believe that the relation between the singular Hele-Shaw problem and the 
Painleve I equation clarifies both entities. It further explains the nature of 
weak solutions of Hele-Shaw problem and shocks in viscous fiuids, and provides 
a transparent, hydrodynamic interpretation of a subtle semi classical description of 
Painleve I, which otherwise is rather formal and complicated. In this paper we give 
a brief account of both sides of the story: the shock-wave solution of the Hele-Shaw 
fiow and the Whitham averaging method for the physical solution of Painleve I, 
even though elements for each of these themes taken separately can be found in the 
literature. We have found that approach to Painleve -equation as an isomonodromic 
deformation developed in Refrs. [ini HZJ [TSl [TSl 1201 [2T] provides the closest connection 
between Painleve -equation to Hele-Shaw hydrodynamics. 

The appearance of the Painleve equation in the Hele-Shaw problem is more than 
a formal analogy. Painleve equation gives an integral Hamiltonian regularization of 
the Hele-Shaw singularities anticipated in Ref. [^. 

Contrary to the Darcy law, the Painleve equation (like a DLA process), has a 
dimensional parameter h. The meaning of this parameter can be understood as 
follows. The physical solution of Painleve I-equation is selected by the asymptote 
u — )■ +\/—T > at large negative times T — )■ — oo. At large positive times T > 0, 
this solution features fast oscillations with a period of the order of h. The period of 
oscillations also changes in time, but slowly, such that the period in units oi h is a. 
regular function at /i — ?■ 0. As such, u{T) does not have a limit at /i — )■ 0, but its 
period (in units of h) does. The Hele-Shaw fiow arises from the Painleve equation 
in the limit when a dimensional parameter h vanishes. It features singularities and 
discontinuities while the Painleve equation a.t h ^ does not. At a non- vanishing 
h shocks are regularized. They appear as spatial narrow regions of width ~ hl^^"^ 
where pressure and velocity exhibit finite, though sharp, gradients. In this sense, 
the Painleve equation provides a regularization of the Hele-Shaw fiow, similar to a 
regularization that provides DLA. 

Emergence of singularities in hydrodynamics is typically caused by an invalid 
approximation made in an original physical problem. Any physical problem of fiuid 
dynamics always carries a dimensional parameter, which regularizes singularities at 
a microscale. As a result, physical processes are smooth, but may exhibit large 
gradients in certain regions of time and space. When the microscale parameter h is 



naively set to zero in the equations, they may become ill-defined and do not have 
physical solutions in the regime when gradients are large. However, if the parameter 
is set to zero in the solution (rather then in the equation), solutions retain their 
meaning but may become discontinuous, exhibiting shocks - weak solutions [9]. This 
is what happens in the Hele-Shaw problem. 

Painleve equation itself is an example of these phenomena. Setting h = in the 
equation reduces it to —2uut = 1, whose solution y/—T makes sense (i.e., remains 
real) only at T < 0. Although at T > the Painleve -function has no formal 
/i — 7- limit, its average over the period of fast oscillations or a motion of minima 
of Painleve -function does (as it is seen on the Fig. \13\). We will show that this 
slow motion is equivalent to the weak solution of the Hele-Shaw problem obtained 
inRefrs. ^^. 

In many important cases, weak solutions show a great degree of universality - a 
large class of regularized problems lead to the same weak solution. It is always 
interesting to determine regularizations that are "minimal" in some sense. In this 
paper, we suggest that the Painleve equation provides the "minimal" regularization 
for generic singularities of the Hele-Shaw problem. It is the only known integrable 
regularization. 

The nature of the physical solution of the Painleve equation suggests a physical 
interpretation of the dimensional parameter h. The solution suggests that, under 
the regularization, the inviscid fiuid consists of particles with a minimal area H. 
Oscillations refiects a discreetness of particles and a "corpuscular" nature of the 
problem. This interpretation links the Painleve I regularization of the Hele-Shaw 
fiow with the DLA process where particles are essentially discrete, and provides an 
analytical basis to study the latter. 

In the next section we briefiy review the Hele-Shaw critical fingers in terms of 
evolution of the real complex curve. It provides a necessary language to link the fiow 
with the Painleve wave-function described in the sequel. Our approach is based on 
isomonodromy property of Painleve -equation. Our approach is somewhat different 
that in the Refrs. p!7l [T9| [20| [21] . In particular we give a new simple derivation of the 



expression for the phase (45) and the next leading oscillatory corrections (63), (64) 
in Sec. refMKB . 



8 



A comment on related studies is in order. In Refrs. [221 123] Hele-Shaw flow has been 
linked to normal random matrix ensembles, and through them to certain asymptotes 
of bi-orthogonal polynomials. Random matrix ensembles provide an equivalent 
regularization scheme of Hele-Shaw flow. Painleve equation can be formally derived 
from this approach ^Hl [21]. In a similar manner Painleve equation appears in 
studies of Hermitian random matrix ensembles |25] and related to their asymptotes 
of orthogonal polynomials [261 ISZ]- From this perspective a singular Hele-Shaw flow 
is linked to a distribution of zeros of orthogonal and bi-orthogonal polynomials and 
eigenvalues of random matrices. In this paper we do not elaborate this approach in 
order to reduce the size of the paper. 

We also comment that an idea that near-cusp structure should be followed by a crack 
propagation into the fluid domain has been suggested in earlier works [2S1 121] • The 
mechanism discussed in these papers seems different from that in [H [2] . 

2. Singularity and viscous shocks 

2.1. Evolving complex curve 

The flow near a singular flnger tip does not depend on details of a boundary away 
from the tip of the flnger. It is sufficient to consider a flnger with a symmetry axis. In 
Cartesian coordinates aligned with a flnger axis (Fig. [I|, a flnger is given by a graph 
±y{x). We recall a description of a viscous flnger in terms of the height function 
Y{X,T) [8\. It is deflned as an analytic function on a flnite part of the complex 
plane X = x + iy cut along the symmetry axis of the flnger, such that the boundary 
values are real and equal to the graph of the flnger Y{X)x=x+±io = ^y{x, T). Darcy 
law can be expressed through the height function [S]: 

StY = -dx(t>. (5) 

where the analytic function 0(X) =■?/' + ip is a complex potential of the flow, V^ is a 
stream function, and p is the pressure. 

The height function Y{X.,T) deflnes a complex curve (or Riemann surface) evolving 
in time. Since at real X = x the height Y = ±y(a;) is real, the curve is real: 

Y{X) = Y(X). 



2.1 Evolving complex curve 9 

In Ref.[8j it lias been shown that in general the complex curve is hyperelliptic, i.e., 
Y^{X) is a polynomial of an odd degree. Its degree is preserved while the coefficients 
evolve in time. 

The flow potential 0(X) has an important interpretation. Darcy law states that 
is a local parameter of the curve. This means the following: the potential (f){X) 
is an analytic function. Locally it can be inverted, so that X, and other physical 
quantities, say, the height function, become functions of 0. Then Darcy law (J5]) states 
that functions {Y{(j)), X{(f))) treated as functions of the potential, are regular and 
single- valued. In the case of critical flow, when {Y, X) is hyperelliptic, the variable (p 
covers the physical plane twice: and —0 correspond to the same X, such that the 
double- valued function Y{X) becomes a single- valued function Y{(j)) on the double 
covering. 

The fact that the potential is a local parameter of the curve is the essence of the Darcy 
law. The law can be formulated entirely in terms of the curve: the curve {Y, X), 
with a local parameter 0, moves in time such that the differential idfl = Y dX + dT 
is closed. In this form Darcy law becomes identical to Whitham equations describing 
slow modulation of fast oscillating non- linear waves p^ [30| [12] . This fact was 
recognized in |15j . 

A consequence is Polubarinova-Kochina equation [3T1 [32] (which in modern literature 
is called "semiclassical string" equation [33] ). Treating Y and X as functions of the 
complex potential 0, the Darcy law ([8| reads 



{X, Y} = 1, (6) 

where the brackets are deflned via {X, Y} = drXd^Y — O^XOtY and time derivative 
is taken at flxed 0. 

The flow is best described by a generating function - a central object of this 
study: 

n{X) = -i I YdX, (7) 



2.2 Integral form of Darcy law 10 

where the integration starts at the tip of the finger e and goes through the fiuid. 
Then the Darcy law reads 

Ot^ = i(j). (8) 



2. 2. Integral form of Darcy law 

Darcy law can be equivalently expressed by integration of the differential dQ along 
various paths in the fluid. Let us integrate ([s]) over a closed path B in the fluid: 

— Im f dn = (f) ^^ X d£ = (f) d^, (9) 

dT J B Jb J b 

-—Re i dVL= - i w-M=- i dp. (10) 

dT Jb Jb Jb 

The imaginary part measures a flux of fluid through the closed path. If the path 
goes around a drain (inflnity), the integral Im ^^ dQ counts a mass of fluid drained 
up to time T. It is QT. In a canonical anti-clockwise orientation this reads 

i-^(i dn = Q>o. (11) 

The real part measures circulation along the path. It follows from the differential 
form of the Darcy law ([I]), valid before a singularity, that the fluid is irrotational 
V X V = and therefore J:Re ^^d^ ~ ^- Therefore Re j^dVL = is a constant. 
Provided that the constant does not depend of the choice of the path, the constant 
can be determined by taking the integral around a drain where § dfl is purely 



imaginary as is in (11 ). Thus Re ^^ dQ = on any contour in the fluid. 



If there are shocks, the height function and dfl jump through a shock. The condition 
Re f^ dQ = for a path crossing a shock must be understood as a sum of the integral 
over a path in the fluid plus an oriented jump of dQ at points where the path crosses 
a shock. A proper way to write this condition is to extend the integration over 
the Riemann surface {Y, X) where the height function and the differential dfl are 
smooth. The physical plane where the hight function suffers a discontinuity consists 
onn patches of different sheets of the curve. Then the above condition reads 

Re (b dQ = 0, any cycle on the curve. (12) 



2.3 (2,3)-singularity and a self-similar elliptic curve 11 



Curves with condition (12) are rather restrictive. We call them Krichever-Boutroux 
curves. They appeared in seemingly different, but related subjects: random matrix 
theory [25], asymptotes of orthogonal polynomials [31] and Painleve transcendents 
[351 ESI EZ]. 



Eqs. (TTp2) combined give a compact and complete formulation of the Hele-Shaw 



problem: find an evolution of a real Krichever-Boutoux curve (condition (11)), of 



a given degree, with respect to its residue (time T) at a marked point (condition 



(ID) 



Contrary to the differential form of the Darcy law ([I]), Eqs. (TTp2) extend beyond 



singularities. They admit a weak solution with shocks. Shocks are curved lines where 
patches of the complex curve covering the physical plane are joined together. At these 



lines, the height function jumps. Condition ( 12 ) being specified for a contour going 
along both sides of a shock reads 

Rediscll|shocks = 0. (13) 

These conditions determine position and motion of shocks. These are the results of 
Refs. [HE]. 

In the next two sections we demonstrate how these conditions determine the most 
generic singularity and also physical solution of the Painleve I equation. 



2.3. (2, 3) -singularity and a self-similar elliptic curve 

Painleve 1 equation is related to the most generic (and the most important) 
singularity - the (2,3) cusp, where locally the graph of the finger is y'^ ~ x^. In 
this case a finger is especially simple [H]: F^ is a polynomial of the third degree - an 
elliptic curve having one branch point at infinity. Fixing a scale and the origin, the 
curve reads: 

-Y^ = AX^ - g^X -g^ = 4(X - e)(X - e2)(X - d), (14) 

where g2 and 5^3 are real time dependent coefficients. One of the branch points of 



the curve ( 14 ) may always be chosen real. We denote it as e. It is the tip of a finger. 



where y{e) = 0. The other two are conjugated ei = 62 (condition (12) excludes a 



possibility of real ei and 62 - in this case disc dfi through the branch cut connecting 



2.4 Smooth flow and a degenerate curve 12 

ei and 62 is real and does not change sign, so that Re f dQ can not vanish). This 
is a real elliptic curve with a negative discriminant (73 ~ 275^1 < 0. The coefficient 



g2 is determined by the drain rate. Eq. (11) gives Q ~ g2- We set the rate Q such 
that 

92 = -12T. (15) 

In this setting the time counts from a cusp event —Y'^ = AX^ occurring at a critical 



time T = 0. A real coefficient g^ is determined by the condition (12). 



This is the simplest curve arising in singular flow. Its distinctive property is self- 



similarity. As it follows from (11),(12) 



Y{X, T) = |T|=^/^r (|T|-i/2x, 1) , (16) 

where Y{X, 1) is a unique universal function with no free parameters. 

The scaling property alone is sufficient to express Q through the height function and 
potential: 

n(X) = -i[ YdX = --(XY -2T(p). (17) 

Je 5 

We briefly describe this curve following Refs. [H [2] . 

2.4- Smooth flow and a degenerate curve 

At T < 0, before the singularity, the flow is smooth, there are no discontinuities. 



This follows also from the condition (12). At (72 > only degenerate elliptic curve 
whose discriminant vanishes, (^f — '^'^dl- The only unknown coefficient is thus 
determined, (73 = — 8(— T)^/^. The two branch points coincide to a real double point: 
ei = 62 = +\/—T > 0. Then e = — 2ei = —2\/—T. The degenerate curve 

Y^ = -4{X~e)(^X+^y , e = -2v^=T, T<0 (18) 

represents a torus with one pinched cycle. The real period of the curve is inflnite, 
the imaginary half-period is J^ °° ^ = vriv 12(— T)~^/^. This is a general situation 
- a continuous flow corresponds to a degenerate curve where all cycles but one are 
pinched [8] (the number of non-degenerate cycles is equal to the connectivity of 
the fluid boundary pTSj). Condition that all cycles located in the bulk of the fluid 



2.4 Smooth flow and a degenerate curve 13 

are pinched foUows from incompressibihty of the fluid. Only in this case the fluid 
potential is analytic. 

Alternatively, we notice two holomorphic functions Y{(j)) and X{(f)) with an 
asymptote Y'^ ~ X^ at large X are polynomials of the third and second degree, 
respectively. The Darcy law in the form of pi) determines the polynomials 

X^-[f)\ein y = 2(^)'-3e(T)(^). (19) 

and yields to a differential equation for time dependence of the coefficients 

edre = -2, (20) 



which prompts (18). Real solution exists only at T < 0. 



2.4-1 ■ Beyond singularity: Genus transition. A finger becomes a cusp when the 
branch point e and the double point — f = ei = 62 merge to a triple point 
e = ei = 62 = 

r2 = _4x3, T = 0. (21) 

After the critical time (72 = — 12T < 0, the curve is not degenerate. The double 
point splits into two branch points Ci 7^ 62 = Ci. No smooth solution is possible. 
The branch points appear in the fiuid as endpoints of shocks. Fig. |4] illustrates the 
genus transition. A pinched cycle of the torus now opens. 



The condition (12) requires that the integral over a new cycle be purely 
imaginary 



/ ^4X3- 



n{ei)= / V4X3 - ^2 a: -g^dX = imaginary. (22) 



Exactly the same equation has been appeared in works of P. Boutroux [35] on 
Painleve I equation. We will refer it as Boutroux condition. This equation can 



be simplified: since Y{ei) = 0, scaling properties (17) yield that at the end of a 
shock ei, pressure 0(ei) = also vanishes. The latter uniquely determines the only 
unknown parameter (73 (T) of the curve. 



2.4 Smooth flow and a degenerate curve 



14 




Figure 3. Stokes- level (dashed-blue) lines Rel^ = and Stokes lines (solid black) 
Im fi = on the fundamental domain of the double covering of the physical plane. 
The blue-shaded region maps to the upper half-plane in the physical plane. The 
other shaded region maps to the lower half-plane. 



Let us parametrize the curve Y{X) by a uniformizing coordinate A as 

X = p{\\ uj{T),u{T)), Y = ip\\\ co{T),u{T)) (23) 

and compute potential of the flow from the Darcy's law and scaling property 



(16) 



- idr^ = (j){X) = 6z (C{X) 



8T' 



(24) 



Here p and C are the elliptic functions of Weierstrafi, whose complex conjugated half- 
periods oj{T) = T~^/^co'(l) and w(T) are yet to be determined. The rhombus-shaped 
fundamental domain is depicted in Fig. [3j 



(25) 



Conditions that pressure vanishes at ei and formula (24) give the deflning 
equation: 

35-3 ^ r] + f] 

where w -|- a) is a real period of the curve and rj = ({u)). 
The solution of this equation is found in 

^3 = \T\'/' 



-7.321762431 



T<0, 
T>0. 



(26) 



2.5 Motion of the finger 15 

uj =T~^/^(0.7426778188 + z 0.6070379047). (27) 

The genus transition is characterized by an abrupt change of g^. The modular 
invariant of the non-degenerate curve 3 l^ 2 = 6.15870507... is a unique 
transcendental number. 

2. 5. Motion of the finger 

The finger is described by the real section of the curve Y,X, i.e., by the graph 



y{x) = Y[X) where X = x < e(T). As it follows from the (18) a finger grows and 
its tips moves to the right as its approaches toward a cusp singularity. After the 
singularity, the finger retreats, moves to the left and triggers two shocks growing to 
the right. The tip of the finger gives the position of the branch cut extended toward 
— CX3. The scaling property shows that tip of the finger moves as \/\T\ but with an 
abrupt change of the coefiicient [H [13] 

,(r) = /-2-^. ^ ^<°. (28) 

^ ' [ -0.553594yT, T > 0. ^ ' 

In the next section we show that the envelope of the Painleve function behaves in a 
similar manner 

«(r)-|-^^(^) = ^^' ^<°' (29) 

^^ 1 e(T) = -0.553594v^, T > 0. ^ ^ 



2.6. Shocks and Stokes level-lines lines. 

Shocks are lines of discontinuity of the generating function f2, the height function 



Y and the potential 0. They are selected by the condition (13), which in the case 
of simple branch points says that shocks are a subset of zero level lines Re i7 = 0. 
We refer to the lines Refi = as Stokes level-lines. There are a total of seven 
Stokes level-lines connected at branch points. They are transcendental, computed 
numerically, and depicted on the right panel in Fig. |4| denoted by Fi, ..., Fy. 

Among the seven Stokes level-lines, only two F3, F2 are shocks. Both sides of the line 
Fi are boundaries of the finger. 



2.6 Shocks and Stokes level-lines lines. 
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Figure 4. Geometry of Stokes-levcl lines (black solid lines) before (left panel) and 
after (right panel) transition. Before the transition four level lines emanate from 
the double point and end at infinity, neither one is a shock. After the transition 
the double point splits into two branch points. Condition Re fi > on the remote 
part of the finger argX = tt, |X| — >■ cxd uniquely determines configuration of shocks 
r3,r2. They are always located in the finite part of the plane. 



Selection of shocks works as follows. We assume that a remote part of the finger 
(argX = vr, |X| — )■ oo) is not affected by the transition. Therefore at large X the 
height function Y{x) jumps only on the finger, i.e., on the negative part of real axis. 
This excludes Stokes level-lines lines r4 5^6,7 ending at infinity. Thus Y jumps only 
on finite Stokes-level lines r2,3. On a shock, Y (as a complex vector) is orthogonal 
to the shock. This uniquely determines the shape of shocks (Fig. |4]). 



The left panel in Fig. |4] represents Stokes level-lines before the transition T < 
There are 9 lines. Four Stokes level-lines r4, Fs, F7, Fg emanate from the doubk 
point. All of them end at infinity. Neither of the Stokes level lines is a shock. 



Having in mind the connection to the Painleve wave-function discussed in the next 
section, it is useful to translate the selection rule in terms of the generating function: 
the signs of ReQ are opposite on both sides of F4^5 6,7, and signs are positive on 
both sides of Fi,F2,F3. This is a general situation: Stokes-level lines Refi = are 
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shocks if Re Q has the same sign (positive) in the proximity of both sides of a shock 

mm 

Ref2|x^shocks > 0, Refi|xeshocks = 0. (30) 

In the next section we show how the same objects appear in the context of non-hnear 
Stokes phenomenon for the Painleve I equation. 



3. Hele-Shaw flow and physical solution of Painleve I 

3.1. Linear problem for Painleve I 

Painleve I equation arises as a result of compatibility of two linear differential 
equations satisfied by the "wave-function" \1/(X, t) - a central object of this 
study. 

Let X and Y, defined through 

X = d^-2u{t), -iY = -2df + 6udt + 3ut (31) 

be two Hermitian differential operators. Then the linear problem reads 

X^ = X^, Y^ = ihdx^. (32) 

The two equations are compatible only iiu{t) obeys Painleve I equation Uu—Gu^ = ht. 
Only then 

[X, Y] = -ih, (33) 



as it follows from (32) 



3.2. Painleve equation as a quantized Hele-Shaw flow 

Equation ^ and (31 ) are to be compared with ^ and ( 19 ) describing a smooth Hele- 
Shaw flow. One identifies the operators X, Y, 0, u with their classical counterparts 
X, Y, (j), — e, where we define := —iHOt = —Qidt. From this point of view Painleve 
equation is seen as a quantization of Hele-Shaw problem, when the brackets in ^ 
are replaced by the commutator ind4|):{,}— )'|[,]. The quantization introduces a 
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new parameter h which regularizes singularities of the flow. It will be sent to zero 
when a solution is found. This point of view has been suggested in [8] . 

As we discussed in the introduction, the limit h -^ exists only before a flow reached 
a singularity. After this point solution oscillates and formally does not have a limit. 
However its average over the period of oscillations does. 

Quantization suggests a regularizing scheme for the Hele-Shaw flow: observables 
in the classical flow are diagonal matrix elements of the corresponding operators 
averaged over the period of oscillations. One may consider the height function, 
flow potential, and velocity of a regularized flow deflned as an expectation value of 
corresponding operators averaged over the period 
{Y) := *-i (thdx) ^, 



(0) := qj-i {-ikdr) ^, (34) 

(v) ■= ^-1 [ikdrdx) ^, 
where O = (period)"^ ^ Odt is an average over the period. These quantities have a 
regular (asymptotic) expansion in yh. 

In the remaining part of the paper we will show that a semiclassical limit of 
the physical Painleve wave-function describes the Hele-Shaw flow before and after 
a singularity. This statement is summarized as: the generating function Vt = 
—i J YdX, (7| of the flow is the Hamilton- Jacobi action of the semiclassical 
approximation 



n = \imh\og^{X,t). (35) 



h^o 



3.3. Physical wave-function and asymptotes at infinity 

The general solution of the Painleve equation is a two-parametric family. Among 
them we must select a physical solution which we corresponds to the Hele-Shaw 



flow by means of (35). A central physical principle we have already used as a 
guideline is: flow around a flnger does not affect the flow at inflnity The height 
function and the potential are holomorphic functions around inflnity (the drain) 



-QiX^/^, Y ~ 2iX3/2 at X ^ oo, lArgX] < vr, as it follows from (19). In 



terms of the wave function, these conditions read: the asymptote of the physical 
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wave-function at infinity is an analytic function in the plane cut along the finger 
ArgX = TT 

/ilog^U^oo = ^U^oo ^ ^X'/^-tX'/\ |Arg|X < TT. (36) 

5 

Also, we must impose two additional conditions. One is reality 



^(X) = ^(X). (37) 



The second is that at t = QT/h — )■ — oo the Boutroux condition (12) holds 



lim ^Re * dlog* = 0, t ^ -oo. (38) 



As a consequence 



uit)-^-^ = V^=J-j>0, t^-oo (39) 

where e(T) is a coordinate of the tip of the finger and e{T)/2 is a position of the 
double point. 

These simple requirements determine both the physical Painleve function and the 
physical Painleve wave-function. 



3.4- Viscous shocks and Stokes phenomena of Painleve transcendent 

The Painleve wave-function is an entire function. It is differentiable everywhere in 
X. However, its semiclassical limit is not. The gradient "^^^hdx"^ is of order one 
in most part of the plane, except some lines where it changes abruptly on a scale 
vanishing with h across these lines. At h ^ sharp gradients become discontinuities 
- seen as shocks in viscous fiow. The origin of these lines is simple: they go through 
the locus of zeroes of the wave-functions. Zeros are separated by h (at the end of the 
line the separation is h"^'^), and at /i — )■ accumulate to a line distribution. Gradients 
through lines of zeroes are large. 

Discontinuity of the semiclassical approximation of analytic functions is known as 
Stokes phenomenon. The locus of zeroes is closely related to Stokes phenomena. 



Shocks occur on some Stokes level-lines selected by the condition (30) and described 
in more details in the sequal. 
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The governing conditions (37) and (36) determine the Stokes multiphers and 
semiclassical wave-function. 

4. The physical solution of Painleve I transcendent 

Among many solutions of Painleve 1 equation 

uu - 6^2 = ht (40) 

called Painleve functions (see e.g., [381 [T7]), we look for the solution which at t — )■ — oo 



corresponds to a smooth Hele-Shaw flow as in (39). This correspondence also 
establishes a relation ^ between the time of the Hele-Shaw flow T and the argument 
of Painleve -function t a relation between times . The correspondence selects a 
one parametric family of Painleve -functions. Following [ll] we call this solution 
"physical" (it is also known as "separatrix" [17] ). 

This solution is scale invariant under the rescaling 

u -^ 6-^/^u, t -^ 5-^'H, ut -^ 6-'^/^Ut, h -^ 5h. (41) 

Below we review some major facts about asymptotes of the physical Painleve - 
functions and derive most of them. 

4-1. Smooth regime t — )• — oo 

The sign of the asymptote u\t^_oo > fixes the coefficients of the asymptotic 
series 

.(t)-^ = -Mp^,,(_t)-5./2, t^-oo, (42) 

fc=i 



where e(T) is given by (39) and t = QT/h. 

The coefficients a^ can be computed recursively, cti = ^- They have a topological 
meaning as intersection numbers of curves on Riemann surfaces [39], and grow as 



«'^~i(i)'^'(fv^)' {2k)\mmm- 
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The asymptotic series alone does not fix the Painleve function. A one-parametric 
family of solutions shares the same series, but differs by the amplitude of an 
exponentially small correction: consider a difference between two different solutions 
sharing the same power series ui and U2- It obeys the equation (mi — U2)tt = 6(mi + 
U2){ui — ^2). Replace Ui +U2 by the power series, and solve the hnear equation with 
respect to Ui—U2- In the first approximation we obtain {ui—U2)tt ~ — 12e(T)(ui— M2). 
The large (negative) time asymptote is controlled by the WKB approximation: 

g-|^-6e(r)(-i) 

Ml -M2 ~ e(T)- 



27rv/-6e(T)(-t) 



This formula is valid in the sector |Argt — 7r| < ^ if (— t) here and in (39) is 



understood as e ^'^t on the complex plane of t cut along the negative real semi- 
axis. 

Along Stokes level lines Argt = vr =1= ^ in the complex plane of t at |t| — )■ 00 
the exponential term becomes of order one. There the asymptote of the physical 
Painleve -function reads [IZl [371 EOl [IE] 



u{t) + \^- — -7^ Z r- Argt = 7rT— (43) 



2n^-6eiT){-t) 

where a real parameter C uniquely characterizes the solution (these formulas are 
written in a form emphasizing slow and fast time dependence e(T) = —2(e^^ht/6Y^'^). 
The amplitude C ± i/2 jumps through the negative part of the real axis. The jump 
is imaginary and is the same for all physical solutions. The jump is linked to the 



large order behavior of the coefficients of the asymptotic power series (42) and can 
be obtained through the Borel transformation technique |18j . 

4-2. Oscillatory regime: t — )■ +00 

As t — ;■ +00 the family of solutions features fast oscillations with a period of order 
h, and has infinite number of double poles clustering at infinity as shown in Fig|2] 
(there are other double poles in the complex plane of time). Semiclassical description 
of the solution at t — )■ 00 is based on a separation of scales between fast oscillations 



4.2 Oscillatory regime: t — )■ +00 
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and the slow-varying period. Separation of scales is controlled by the parameter h, 
or large T. 

At positive large time, the physical Painleve functions are asymptotically Weierstrass 
elliptic function 



u{t) -p{^t- eiT)\u{T),co{T)] , T = h^+oo. 



(44) 



The complex conjugated periods of the Weierstarass function and the phase 9{T) = 
slowly depend on time and in the leading approximation do not depend on h 



The periods are determined by the Boutroux condition [35] (see ([49|) below) which 



happens to be identical to the integral formulation of the Darcy law (12), (22) and 



is given by (27). Eq.(44) valid at complex time in a sector bounded by Stokes lines 
Argt 



±1- 



A value of the phase depends on the solution chosen. If one choses to characterize 



the Painleve -function by the constant C in (43), one can also characterize the phase 



by C. Relation between C and 9 links asymptotes of Painleve -function at T < and 
T > 0. It is called connection formula and obtained in Refrs. [371 [IS]- Connection 
formula could be best expressed through the monodromy data $ - a (locally) time 



independent constant defined below in Sec. 5.4 
ne = ^LU + $0;, C = -Im e"^* 



2Ree 



j<E> 



-1. 



(45) 



The phase 6 is understood through initial data for Painleve equation. As it follows 



from (44), at times 
4 



9 + {2n + l){u + u!), 



(46) 



where n is a large positive integer, the Painleve -function reaches the minima 
Ut{tn) = as in Fig. |2} At those points the value of the solution u{tn) -^ e(T), 
where e(T) is a real branch point of the Krichever-Boutrox elliptic curve ([28|. 



Scale invariance (41) implies that the periods behave as \T\ ^/^ (27). Taking this 



4.3 Adiabatic invariant and Boutroux condition 



23 



into account the formula (44) reads 



u{t) 



6 



P 



4 fhtV^ 

5M-J -^(1) c.(l),a;(l) 



t -^ +CX), (47) 



where 6'(1), oj{1) are just numerical constants. 

We will derive all the formulas of the oscillatory regime below. Connection formula 



(45) requires a more subtle analysis and we do not discuss here. 



4-3. Adiabatic invariant and Boutroux condition 

Painleve equation can be interpreted as mechanics of a particle with a coordinate 
u moving in a potential —2u^, under slowly changing force ht. At t — ;■ — 00 the 
particle with a vanishing velocity is located in the unstable position where the force 



equilibrates the potential as it follows from (39). From there as time progresses the 



force no longer equilibrates the particle. The particle moves to infinity, but returns 
from there at a finite time, performing an oscillatory motion. 



In oscillatory regime the r.h.s. of the Painleve - equation (40) is a "slow time" T 



barely changing over many periods of oscillations. Treating the r.h.s. of (40) as 



a parameter, the equation becomes an identity for the Weierstrass elliptic function 
p" — 6p^ ~ ~^' where (72 = — 12T. This gives (47), where the Weierstrass constant 
f73 and the phase 6 remain undetermined. 

Adiabatic nature of modulation of oscillations helps to determine the Weierstrass 
constant g^ and the periods of oscillations. 



The energy of the oscillatory motion is 

H = P^ -Au^ + g2U, 
where P 



(48) 



Ut 



y^Au'^ — g2U + H is the momentum of oscillations. 



The Weierstarass parameter g^ is a value of the energy (48) evaluated on (46): 
H — i- —(73. The adiabatic invariance of the motion f Pdu does not depends on 
time. At earlier time when the oscillations do not yet begin - the adiabatic invariant 
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is null. It remains null at all time. 



® a/4m3 



g2U - g^du = 0. 



(49) 



The integration in this formula goes over the period of oscillations. 



The fact that the adiabatic invariant vanishes also means that at early time the mean 
value of the energy (48| Hy^_^ = —g^ — >■ (5^2/6) = e^ = — 8T^/^ is such that the 
curve iu^ —g2U — g^ = 4:{u — e){u+ |)^ degenerates as the period of oscillations tends 
to infinity. 



Eq.(49) is called Boutroux condition [35]. It determines the mean energy H = —g^ 



oscillations and the Weierstarass parameter (73 and subsequently determines the 



periods. Boutroux condition (49) is equivalent to the condition (22) obtained from 



the weak solution to the Hele-Shaw problem. The calculations and the results 



(26), (27) were described in Sec. 2.4.1 



Direct integration of the Painleve equation also yields the Boutroux condition. First 
we write the Painleve - equation in the form which emphasizes a non-conservative 
nature of motion 



Ht = -2hu, 



(50) 



where H is a. time-dependent Hamiltonian of the mechanical system is in (48). Then 
compute and compare the means (average over oscillations) of both sides of the 
equation. In the leading approximation g^iT) = g^{l)T^^'^ depends only on "slow 
time". Theieioie OtH ^ -^ = 12|^. The integral over the period is 2-12(a;-Fw''^ 



'292' 



It must be equal to the integral over the period of the r.h.s. of (|50|). With the help 

1 pt+2{uj+u)) 
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of (47) we obtain u 



Ul-\-U) 



/. 



u{t)dt 



v+v 

LJ+Gj ' 



This yields Eq. (25) which is 



equivalent to (49) 



In the next section we refine these results through the monodromy of the associated 
linear problem. In that way the Boutroux condition emerges exactly in the form of 
the integrated Darcy law. 
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5. Semiclassical analysis of Painleve I transcendent 

5.1. Linear equation for the Painleve -wave-function 

A starting point for the semiclassical analysis of the wave-function is a commonly 
used realization of the linear problem [H], equivalent to Q or (32). Differentiating 



the first equation (32) by time, one expresses d^^ in the second equation through ^ 
and dt^ . Then differentiating the result by time again, one expresses d^^ through 
^ . The result can be written as a linear matrix differential equation for the column 
vector "^ = (\l/, \l/t) treated as a function of space X 

ihdx-^ = Y£, (51) 

2x2-matrixY explicitly depends on u, Ut and t (we express Uu through the Painleve 
equation) 

-iY = utas + 2 (X - m) a+ - 2 ((X + 2u){X - u) + 3{u^ + T)) a_, 
where as, a± are Pauli matrices. 



This linear equation alone, without a reference to (32), already determines that u{t) 
is a Painleve function, if one demands that monodromy data does not change in time, 
while the matrix Y does. This is the isomonodromic deformation interpretation of 
the Painleve equation [TCI EI] • 



The linear equation (51 ) has one irregular singular point at infinity, and three turning 



points where det Y = 0. Solutions are entire functions. 



5.2. Modulated spectral curve 

The eigenvalues ±Y?i of the traceless matrix Y are two branches of a spectral elliptic 
curve 

- Y^ = -detY = 4X^ - g2X + H. (52) 



Here H is given by (48), and, again, g2 = — 12T = —2ht. Another form of the curve 
is also useful 



-F;f = 4(X + 2v/=T) (X - V^y -4(u + 2^^^ (u - v^=t)' + 



2 
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It is organized sucli tliat at T — )■ — oo one easily recovers a degenerate curve of the 



Hele-Shaw flow (18). The last two terms vanish in that limit. 



We emphasize that Y^ features fast oscillations with oscillations of energy H. The 
average over the period of the oscillations H ^ H = —g^ revials a slowly evolving 
modulated spectral curve Y{X). 

5.3. Semiclassical wave-functions 



The spectral curve determines the leading order of the semiclassical approximation of 
the wave-function, while eigenvectors of the matrix Y determine the sub-leading term 
[§] [12]. Two independent basis WKB-functions correspond to each branch of Y 



^^^^(X,t) 



± 



X 



u 



Y 



1/2 



2 J. 



g 2 JX (X'-u)(iY(X'))Q^-K 



±iQfi 



(53) 



where rffi/j = —iYndX. The formula holds away of a \/h vicinity of three turning 
points where Y^ vanishes. The function ip^ matches the asymptote of the wave- 



function (36) at X — )■ cxD, which is computed directly from (51) 

H/8h 



vl/(X,t) = ^. 



WKB 



X 



+ 



i[lx2-txh 



(54) 



In further approximation we replace the oscillating H by its mean H and call it —gs. 



We already know that g^ is determined either by (25), or (49), but we do not assume 
it at this point. 



Semiclassical form of the Painleve -function can be obtained from (53) under 



requirement that the WKB basis is locally single valued functions on the double- 
covering of the curve (the rhombus in Fig. [3]) with branch points at inflnity and at 
turning points ei,2- 



We set X = p{X), —iY = p'{X) uniformmizing the curve (52) and choose A to be in 



the shaded domain of the rhombus in Fig[3J Let us also set u{t) to be the Weierstrass 
function u{t) = p(r) with the same periods as X and Y, where r is some function 



§ The leading term of the WKB sohitions is exp J (^^^ri ^ {U ^dU)ii), where U is a. matrix 
diagonaUzing Y : U^^YU = Yna^. Each solution corresponds to a branch of Yn- 
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of time. Then, with a help of the identity — (a) (r) ~ ~Ci^ ~ ^) + Cl-^ + ''")" 2C(r) 
we compute -f /^(x^zSw) = f /a° Ra^ = 2#M (log ^ " 2AC(r: 
Another identity reads X — u = p{X) — p{t) = '^2(x)a2(T) ■ Combining the pieces 
we observe that the WKB -wave-functions are meromorphic only if Ut = p'{t). The 



latter determines the argument of p function as in (47) 

4 

r = -t-e. 

leaving the phase ^ as a parameter of solution. 
Summing up, the WKB solution reads 

<™(A.<) = ^^.->'^f^V AeB, (65) 

a(A)c7(r) \VY J 

where A is chosen to be in the shaded domain B of the rhombus Fig. [3j Here a and 
( are elliptic functions of Weiertstrafi. We grouped terms such as to emphasize the 
slow time dependence of the spectral curve Y{X) and fast oscillations due to "fast 
time" dependence of eigenvectors of Y. 

The second component of the WKB basis can be obtained as a result of rotation of 



( 55 ) around infinity [U] 



^M^^^(A, t) = iV^^^^(e^-A, t), XeB, (56) 






^wKB^X^2m^ = e^V-- ''''(^) (57) 

At moments t„ in (|46|, the WKB wave-functions are especially simple 



^^^^'^■'")-^'/ (x4)^-e/ ^'" (««' 



5.4- Monodromy data 



The basis of WKB solution is quasi-periodic 

^WKB^^ + 2io) = -^^^^(A)eT^* 
^WKB^^ + 2io) = -^^^^(A)e 

Complex Bloch factor $ constitutes the monodromy data 



WKB{\ I n,-,\ — \T,WKB{\\„±i<!> (^9) 
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The fundamental fact about the monodromy data is that they do not depend on 
time if u{t) is a Painleve -function. The reverse is also true - if $ does not depend 
on time, then u{t) is a Painleve -function. This statement in full generality was 
proven in Ref . [16] . For the case of Painleve I equation the time independence of the 
monodromy data can be seen from the following simple argument: Wronskians of 
Q computed out of pairs of solutions which WKB asymptotes in a given sector are 
^7 ^^(A), ^^^^(A) and Vf ^^(A), ^5^^(A+2a;) are related by the Bloch multiplier. 
Both do not depend on time. In a scaling solution Bloch multipliers also do not 
depend on h - they are two complex conjugated numbers which could be chosen to 
characterize the Painleve -function. As it has been shown in [ITj physical solution is 
selected by the condition 2Ree** = —1. 



A jump of WKB wave-functions (59 ) on branch cuts can be expressed through Bloch 
factors. With a help of (g we write ip^^^iX + 2u) = iip^^^^-X - 2u) = -e"^*. 
The argument —A — 2ci; is a rotation of a point A in the rhombus by vr around —Co. In 
the X-space it represents a rotation by 2tt around the turning point 62- Thus 

^f ^^(ei + Xe2-*) = ie^'^^p^^^{el + X) 
^f ^^(62 + Xe-2-^) = ^eT*V^^^^(e2 + X), (60) 



Together with (56) these equation describe jumps of the WKB-basis on branch 
cuts. The factors Si = — ie"** and s_i = — ze** are often called Stokes multipliers 

5.5. Monodromy data Krichever-Boutroux condition 



An explicit form of the wave-functions (55) contains all essential information about 
Painleve -function. We have 

- z<l = 2 (a;C(r) - m) - t (^fi(A + 2a;) - fi,(A)) (61) 

h 



2uj{ar)-~]+2e^, 



where we used 



fift(A + 2uj) - fift(A) = ^92 Ui + ^uj ) (62) 
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and standard formulas a{\ + 2lo) = -e2''(^+'^)(T(A), C(A + 2cu) = C(A) + 2r/. 



Taking a real part of ( plj ) we obtain Re ^ dVLj^ = Re (fi^ {\ + 2u + 2a;) — ^^(A)) = 
2h{uj + id) (C(t) — zj^^) ~ ^Ini^. Averaging over one time period cancels the 
oscillatory term in the r.h.s.. Hence 



Re ® d£tn 



-/ilm$. 



(63) 



The limit /i — ;■ yields the Krichever-Boutroux curves (12). 



Another form of this result is a correction to the mean energy: 

hh Im<l> 



H 



-93 



3 uj 



u 



(64) 



It follows from averaging the real part of (61) over the time-period 



Finally, the relation between the phase and the monodromy data (45) vr^ = $a; + $a; 



follows from (61 ) - multiply (61 ) by cj, take the imaginary part and use the Legendre 



relation rju — rjuj 



■TTl 

2 ■ 



One checks directly that in the leading order (61) is consistent with (50) only if the 



monodromy data are time independent. Differentiate (61) by the fast time and take 



into account that all terms in (61) which do not scale with a slow time. Eq.(50) 
averaged follows. 



6. Semiclassical wave- function 



Now we are ready to discuss the WKB-approximation of the physical Painleve -I 



wave-function determined by the condition (36). The WKB form depends on the 
sector 5* bounded by two consecutive Stokes lines Imr2(X) = (blue-dashed lines 
on Fig. |5| and a real negative axis. 



WKB, 



^(X) = A+{S) i/j'^'^'^iX) + A^S) ^!^^^(X), X eS^ 



(65) 



Away from Stokes level-lines where Re Q,{X) = 0, only one term in (65 ) is relevant. It 
is the dominant term: either ijj^^^ if Re ^2 > 0, or ^^-^^ if Re fi < 0, provided that 
its coefficient does not vanish. If the coefficient of the dominant term vanishes, 
the WKB approximation is given by the recessive term ipY^^ if Ref2 < 0, or 
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^WKB if j{^gf2 > 0. On Stokes-level lines, both terms in (65) oscillate and are 
comparable. 



If both coefficients in (65) are non-zero in some sector, the physical wave-function 



has a set of zeros located on a Stokes-level line bisected this sector. Zeros accumulate 
toward turning points - end points of the lines. In Hele-Shaw hydrodynamics (Sec. 
[2]), some of these particular Stokes-level lines are located in the fluid. They appear 

as shocks. 

In order to complete the semiclassical analysis we have to find the Stokes coefficients 
A±[S) for the physical Painleve -functions. They are expressed through Bloch 
multipliers $ and depend neither on time nor on h. 




Figure 5. Stokes lines Imf2 = (blue-dashed lines) before (left panel) and 
after (right panel) transition. Stokes -level lines are black solid lines. Before the 
transition four Stokes lines emanate from the double point. Stokes phenomena 
occurs only on Fi representing the fluid boundary. After the transition the double 
point splits into two branch points, each emanating three Stokes lines. Stokes 
phenomena occur on shocks r3,r2 where the Painleve -wave- function accumulates 
zeros. 
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6.1. Stokes coefficients 

In further analysis, Stokes level-lines (also called anti-Stokes lines), Stokes lines play 
an important role. They have been already discussed in Secj56| We recall the major 
facts. 

Stokes level-lines are curved lines determined by the condition Rer2(X) = 0. On 
these lines the basis WKB-functions oscillate. There are seven level lines as indicated 
on the right panel in Fig. Isj Three of them (ri_2,3) are distinct: YleVt approaching 
to these lines stays positive on both sides. On other lines Re f2, changes sign. Being 
positive in the sector bounded by r7^3^2,5), Refi becomes negative in adjacent sectors 
bounded by r5^4 and in the sector bounded by Fgj, and positive again in the sectors 
bounded by r4 2,i) and Fi ag). 

Stokes lines are the lines where IvuVt^X) = const, and connected to a branch point. 
There are nine of them as shown in the right panel of Fig. |5j All Stokes lines end at 
infinity. 

The semiclassical wave- function is analytic between adjacent Stokes lines separated 
by a Stokes level line, but may be discontinuous passing through a Stokes line. 



Condition (36) requires that at large X the semiclassical approximation suffers a 
discontinuity only at a fluid boundary ArgX = tt where the wave-function rises as 
|X| — > oo. This means that on the both sides of the fluid boundary the wave-function 
possesses a component ip^^^ which happens to be dominant in the sectors bounded 
by Fi,F4 and Fi,F6. In other words A^ 7^ in these sectors. We set it to be 
1. Passing the Stokes-level line F4 the component 1^+ becomes recessive. Therefore 
A^ = to the right of F4, and thus holds everywhere in the sector bounded by S'^ 
and S'ni'. The symmetry with respect to complex conjugation insures that A_ = 
also in the sector bounded by Siv and Sj". Thus 

\[/ ^ ^^^^ between Siy — Si" and Sy — S'ni'. (66) 

In the sectors adjacent to Fi,F2,F3 the WKB approximation of ijj^^^ is dominant. 
There is no information so far about A^ there. Monodromy property of WKB 
functions determines this coefficient. 

Passing through the Stokes-level Fi, the term A^ip^^^ becomes dominant. With 
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a help of (56) we evaluate it at a remote point Xi = |Xi|e*^'^"*"'''*, |Xi| — )■ oo on the 



lower side of the Fi {6 > 0). We have —iA_ilj+{i\Xi\e''^ 7r+<5))_ g^^ ^j^g coefficient of 
the dominant term is always 1. Thus A__ = i in a sector between Fi and Sm. 

Similar arguments determine the Stokes multipliers between Stokes lines Sm, Siw 
and Si, Sif. Summing up 

' ^P^^^ + z^5^^, between Fi, ^i„, 
_ I ^WKB _ ,^r^^WKB^ between Snv, S^i, ..„. 

^ V^r^'^' + ^e^*^!^^'', between^!,, ^i, ^ ' 

^WKB ^ between Sy, Smi. 

The WKB-approximation of the physical wave function in sectors in the lower half 
plane is achieved by complex conjugation. A special feature of this solution is that 
at infinity, everywhere in a fluid (a complex plane cut by Fi) \l/ k. ip^^^ and Eq. 



(36) holds. 



Situation is different at T < 0. Configuration of Stokes lines and Stokes-level lines is 
depicted in the left panel of Fig. |5| In this case the only Stokes-level line where ^Q 
stays positive on both sides is Fi. Similar analysis yields a very different result 

^ ^ . ^f ^^ + i^P'^'"'', between Fi, Suv, .gg^ 

ip^^^, elsewhere in the upper half plane. 



7. Semiclassical wave-function and the Hele-Shaw flow 

7.1. Darcy Law as an average over the period 

The Darcy law follows from the semiclassical wave-function. As it has been described 



in Sec. 3.2|one recovers the hydrodynamics by averaging oscillating quantities over 



the period. Let us see how does the wave-function change over the time period 

2{UJ + L0). 

Away from sectors adjacent to F2, F3 and a fluid boundary Fi the WKB wave- 
function is equal to ip^^^. Let us compute an average /ilog\E' over fast 
oscillations. 



7.2 Locus of zeroes of the physical wave function and shocks 33 



First we notice that the term placed in (|55|) outside of the brackets J^A^n / A e ''^'^•'^ is 



time-periodic. The only change over the period comes from the slowly varying part 
y^^/^es^ in (p5J). Thus the leading contribution into hlog"^ is il. Eqs. (34) recover 
the weak solution of the Hele-Shaw problem. 

7.2. Locus of zeroes of the physical wave function and shocks 

Physical wave-function oscillates on Stokes level lines. Three Stokes level lines ri.2,3 
are special. There, the oscillating WKB wave function possesses two oscillatory 



terms (67), and, therefore, has zeroes. Position of zeroes Xk on these level lines is 



determined by Bohr-Sommerfeld conditions 

f thnik - I), Xk e Pi, 

nn{Xk) = l thikn + l+'^^), X, eP2, (69) 

[ ihikn + l-l^), X, GP3, 

where we evaluate Qn from the above of the cuts. We observe a simple pattern 
selecting the locus of zeros. These are the lines on which Re Q vanishes being positive 



at both sides as in (30) 



Crossing these lines the continuous wave-function experiences an abrupt change on a 
scale vanishing at /i — )■ 0. In that limit the wave-function suffers a discontinuity. The 



height function, potential, etc., (34) obtained by differentiation of the wave-function 
follow suit. In the hydrodynamic limit /i — !■ the lines of accumulating zeros in 
the fluid appear as viscous shocks, each zero appear as micro-vortex as has been 



discussed in (Sec 2.4.1) 



8. Discussion and conclusions 

We have demonstrated that viscous shocks in singular Hele-Shaw flow are intimately 
related to non-linear Stokes phenomena of Painleve I - equation. Stokes-level lines in 
a complex plane accumulating zeros of the Painleve -wave- functions evolve as viscous 
shocks. 

Hydrodynamics quantities have a transparent interpretation in terms of evolution of 
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semiclassical approximation of the physical Painleve - function. In particular, each 
zero of the holomorphic Painleve -wave-function is understood as a vortex in the 
fluid. 

Semiclassical approximation of Painleve -wave-function at T > features fast time 
oscillations and also fast spatial oscillations along shocks. Hydrodynamics can be 
read from these objects either through averaging them over periods of oscillations 



(34), or, alternatively, by considering discrete moments of time (46) corresponding 



to half-period where Ut = 0. This sequence of times depends on monodromy data, 
while hydrodynamics does not 

4/5 



(if)"'' ((^ + Re<l>)Reu;(l)-Im<l>Imu;(l 



'5ft\4/5 
>12/ 

^ (fl ■ 0.7427)^^'' (n + Re $ - Im* ■ 0.8173)' 
Here n is large positive integer. At these moments the wave-function reads 

X - e(T) 



,WKB/ 



|*''^"^(X)|'~ 



X2 + e(T)X + 
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4e(T) 



g|.Im(Xy-2T<^{X)) 



(70) 



where 



Y{X) = v^-4X3 + 12TX 



4>{X) = 6 



Je{T) ^ 



93, 

dX' 



TJ Y{X') 
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The numerical value of g^ and e(T) are given in (26) and (28) respectively by 
gs ^ -7.321762 T^/^ e(T) ^ -0.553594yT. 



Eqs. (34) yields hydrodynamics results for potential and velocity of the flow 



supplemented by semiclassical corrections. 

Discrete nature of time and discrete structure of shocks supports a "corpuscular" 
interpretation of singular flow. One may associate a space between zeros in a shock 
line with a particle aggregating on a shock. A particle arrives to a shock every time 
interval equal to h"^^^. 

Finally we comment that the present analysis does not discriminates within a one- 
parametric family of physical solutions characterized by the Bloch factor $. We 



35 



conjecture that a unique solution will be selected by the requirement that Painleve - 
function is associated with some stochastic process. This requirement may select 
$. An interesting candidate is Im$ = 0. Then the Bloch factor is a pure phase. 



According to (45) $ = 7r/6 and |C| = a/3/2. 
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